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We study the dynamics of a nanomechanical resonator (NMR) subject to a measurement by a low- 
transparency quantum point contact (QPC) or tunnel junction in the non-Markovian domain. We 
derive the non-Markovian number-resolved (conditional) and unconditional master equations valid 
to second order in the tunneling Hamiltonian without making the rotating-wave approximation and 
the Markovian approximation, generally made for systems in quantum optics. Our non-Markovian 
master equation reduces, in appropriate limits, to various Markovian versions of master equations 
in the literature. We find considerable difference in dynamics between the non-Markovian cases 
and its Markovian counterparts. We also calculate the time-dependent transport current through 
the QPC which contains information about the measured NMR system. We find an extra transient 
current term proportional to the expectation value of the symmetrized product of the position and 
momentum operators of the NMR. This extra current term, with a coefficient coming from the 
combination of the imaginary parts of the QPC reservoir correlation functions, has a substantial 
contribution to the total transient current in the non-Markovian case, but was generally ignored in 
the studies of the same problem in the literature. Considering the contribution of this extra term, 
we show that a significantly qualitative and quantitative difference in the total transient current 
between the non-Markovian and the Markovian wide-band-limit cases can be observed. Thus, it 
may serve as a witness or signature of the non-Markovian features in the coupled NMR-QPC system. 

PACS numbers: 03.65.Ca, 03.65.Yz, 42.50.Lc 



I. INTRODUCTION 



Recent advances in nanotcchnology have enabled the 
fabrication of very small quantum electronic devices 
that incorporate mechanical degrees of freedom, called 
nanomechanical systems 1 " - —. The interplay between elec- 
tronic and mechanical degrees of freedom has gener- 
ated interesting dynamical effects^—. These advances 
have also opened a new avenue to technology of high 
precision displacement measurement using electronic de- 
vices, such as quantum dots, single electron transistors 
(SET's), or quantum point contacts (QPC's) 1 ^— . Ex- 
periments using SET's and QPC's have demonstrated 
displacement detections of a nanomechanical resonator 
(NMR) with sensitivities close to the standard quan- 
tum limiki2r— . Similar problems of a two-level system 
measured by QPC's or SET's have also attracted much 
attention 2 ^— theoretically and experimentally. 

The transport properties in the nanostructure elec- 
tronic devices are often studied theoretically in the 
wide-band limit (WBL) and under the Markovian 
approximatio n 12 ' 17 " —. The WBL approximation ne- 
glects an important fact that electron tunneling ampli- 
tudes and also electrodes' densities of states are in gen- 
eral energy-dependent. The Markovian approximation 
assumes that the correlation time of the electrons in the 
electrodes (reservoirs) is much shorter than the typical 
system response time. These approximations may not 
be always true in realistic nanostructure devices. Hence, 
a recent development in quantum nanostructure elec- 
tronic transport has been devoted to the study of the 



non-Markovian effects on the electron transport proper- 
ties, taking into account the energy- dependent spectral 
density of electrode s 33 ' 35 " — . In this paper, we investi- 
gate the dynamics of a NMR subject to a measurement 
by a low-transparency QPC or a tunnel junction in the 
non-Markovian domain. This problem has been exten- 
sively studied in the literature under various conditions 
and approximations 1 ^ - — . In Ref. [Til a master equation 
of the reduced density matrix of a NMR was obtained 
for zero-temperature QPC reservoirs (electrodes) in the 
high-bias limit. The master equation presented in Ref.[2(il 
included not only the effect of the QPC reservoirs in the 
high-bias limit but also the effect of a high-temperature 
thermal environment on the NMR. The master equa- 
tion derived in Ref. [U was claimed to be applicable for 
a broad range of QPC temperatures and bias voltages 
and for arbitrary NMR frequencies. However, the results 
presented in these papers^— were under the Marko- 
vian approximation and without considering the energy- 
dependent spectral density of electrodes. In this paper, 
we take these into account and derive a time-local (time- 
convolutionless) non-Markovian master equatio n 43 ' 45 " — 
that reduces, in appropriate limits, to various Markovian 
versions of the master equations in these papers^—. 
We find considerable differences in dynamics between 
the non-Markovian case and its Markovian counterpart 
in some parameter regime. We also calculate the time- 
dependent transport current through the QPC which 
contains information about the measured NMR system. 
We find an extra transient current term proportional to 
the expectation value of the symmetrized product of the 
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position and momentum operators of the NMR. This ex- 
tra current term, with a coefficient coming from the com- 
bination of the imaginary parts of the QPC reservoir cor- 
relation functions, has a substantial contribution to the 
total transient current in the non-Markovian case and 
differs qualitatively and quantitatively from its Marko- 
vian WBL counterpart. But it was generally ignored in 
the studies of the same problem in the literature^—. 
Considering the contribution of this extra term, we show 
in this paper that a significantly qualitative and quan- 
titative difference in the total transient current between 
the non-Markovian and the Markovian WBL cases can be 
observed. Thus, it may serve as a witness or signature 
of the non-Markovian features in the coupled NMR-QPC 
system. 

The paper is organized as follows. In Sec. [TT1 
we describe our NMR-QPC model. In Sec. Hill 
we derive a time-local (time-convolutionless) non- 
Markovia n. ' — number-resolved or n- resolved 
(conditional) 20 ' 21 ' 24-27 ' 29 ' 31 ' 32 master equation of the 
density matrix of the NMR subject to a measurement 
of a QPC detector and an influence of a thermal bath. 
Our non-Markovian equation is valid for arbitrary bath 
temperatures, electrode reservoir temperatures, bias 
voltages and NMR frequencies as long as the approxi- 
mation used in our approach, namely the second-order 
perturbation in the system-bath and system-reservoir 
coupling strengths, holds. In Sec. IIV1 we present the 
unconditional non-Markovian master equation and show 
that the unconditional non-Markovian master equation 
we obtain reduces, in appropriate limits, to various 
Markovian versions of the master equations in the 
literature. In Sec. [V] we calculate the non-Markovian 
expectation values of the NMR dynamical variables. 
The time-dependent transport current through the QPC 
which contains information about the measured NMR 
system is calculated in Sec. ED We follow Ref. H]] to 
categorize the non-Markovian average current into sev- 
eral physically distinct contributions. We find an extra 
transient current term that has a substantial contribu- 
tion to the total transient current in the Non-Markovian 
case. Numerical results together with discussions are 
presented in Sec. I VIII A conclusion is given in Sec lVIIIl 

II. HAMILTONIAN OF THE NMR-QPC MODEL 

In this section, we describe the model of a NMR that 
is subject to a measurement by a low-transparency QPC 
or electric tunnel junction^ - — and is under an influence 
of a thermal environment (see Fig. [1} . In this model, the 
NMR is considered as a quantum harmonic oscillator, 
and the thermal environment and the QPC electrodes 
(leads) are treated as an equilibrium bosonic bath and 
non-equilibrium fermionic reservoirs, respectively. By 
considering the NMR linearly coupled to the QPC, the 
Hamiltonian can then be written as 

H = H S + H B + Hi, (1) 




QPC 

FIG. 1. (Color online) Schematic diagram of a nanomechani- 
cal resonator (NMR) coupled to a thermal reservoir and mea- 
sured by a quantum point contact (QPC) detector. 

where 

Hs = ^ + (2) 
H B = H Xcads + ]T + \m n wlqi) (3) 

with 

ffloads^ #lead= J2 J2 e k4,k C l,k ( 4 ) 

l=S,D l=S,D k 

and 

Hj -^tunneling + ^ Kiq n x, (5) 
n 

with 

^tunneling = ^(?fc<f + XkqX)Cs,k C D,Q + H.C.. (6) 
k . q 

Here Hs represents the Hamiltonian of the NMR system, 
and m and ui are the mass and the (rcnormalizcd) nat- 
ural frequency of the NMR, respectively. Hb represents 
the Hamiltonian for the left and right leads (reservoirs) 
of the QPC and the thermal (bosonic) bath. The thermal 
bath in Eq. ([3]) consists of a large number of harmonic os- 
cillators with masses m n and frequencies w„ , respectively. 
In Eq. (d]), c^k and e l k are, respectively, the fermion (elec- 
tron) reservoir annihilation operators and energies with 
wave vector k for the left (source) or right (drain) lead 
of the QPC. The interaction Hamiltonian H], Eq. (O, 
contains two parts: the first term describes the tunnel- 
ing Hamiltonian of the electrons through the QPC and 
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the second term describes the interaction between the 
NMR and the thermal environment. In Eq. ([6]), the bare 
tunneling amplitude between respective states with wave 
vectors k and q in the left and right leads (reservoirs) of 
the QPC is given by Tk q , and H.c. stands for the Hermi- 
tian conjugate of the previous term. So the interaction 
between the NMR and QPC introduces an effective tun- 
neling amplitude^ from T^ q — > T^ q + XkqX in Eq. ([5]). 
The NMR and each of the thermal bath oscillators inter- 
act bilinearly through their respective position operators 
as shown in the last term of Eq. ([5]), and A n is the cou- 
pling strength. 



III. NUMBER-RESOLVED QUANTUM 
MASTER EQUATION 

Non-Markovian dynamics usually means that the cur- 
rent time evolution of the system state depends on its 
history, and the memory effects typically enter through 
integrals over the past state history. However, the non- 
Markovian system dynamics of some class of open quan- 
tum system models may be summed up and expressed as 
a time-local, convolutionless form 5 ^ where the dynamics 
is determined by the system state at the current time 
t only. This time-local, convolutionless class of open 
quantum systems may be treated exactly without any 
approximation. The quantum Brownian motion model 
or the damped harmonic oscillator bilinearly coupled to 
a bosonic bath of harmonic oscillators 5 ^— is a famous 
example of this class. The pure-dephasing spin-boson 
model 5 ^— also belongs to this class. The non-Markovian 
effects in the master equations are taken into account by 
the time-dependent decoherence, damping and/or diffu- 
sion coefficients instead of convolution memory integrals. 

The perturbative non-Markovian open quantum 
system theory may also be categorized into two 
classes: the time-nonlocal and the time-local (or time- 
convolutionless) methods. This has been discussed ex- 
tensively in the literature 4 ^ 48 ' 51 . In the quantum mas- 
ter equation approach, after the Born approximation, the 
master equation of the reduced density matrix p{t) of the 
system could be an integro-differential equation and thus 
nonlocal in time. In the interaction picture, the master 
equation in this case can be written a a 48 ' 65 ' 66 



satisfies the condition of 



dp(t) 
dt 



1 



:Tr f 



dt' 



H T {t), ff z (t'),p(«0«>-Ro 



(7) 



where p{t) is the reduced density matrix of the system 
and Hj(t) is the system- reservoir interaction Hamilto- 
nian in the interaction picture. In obtaining Eq. ([7]), the 
assumption of the initial total density matrix in the un- 
corrected (factorized) form of pr(0) = p(0) <8> Ro (with 
Rq being the reservoir density matrix) and the assump- 
tion that the system-reservoir interaction Hamiltonian 



Tr R [H!(t)R Q } = 0, 



(8) 



to eliminate the first-order term in Hi are made. How- 
ever, it can also be shown that another systemati- 
cally perturbative non-Markovian master equation that 
is local in tim o 43 ' 45 ~ — can be derived from the time- 
convolutionless projection operator formalism^ 5 . - — or 
from the iteration expansion method 49 . Under the simi- 
lar assumptions of the factorized initial system-reservoir 
density matrix state and Eq. ([5]), the second-order time- 
convolutionless master equation can be obtained as 



dp(t) 
dt 



1 



Tr, 



dt' 



H r {t), #i(0,p(t)®iZo 



(9) 



We note here that obtaining the time-convolutionless 
non-Markovian master equation perturbatively up to 
only second order in the interaction Hamiltonian is equiv- 
alent to replacing p(t') with p(t) in Eq. $jp!L2L2±.. One 
may be tempted to think that the second-order time- 
nonlocal master equation ([7]) is more accurate than 
the second-order time-local (time-convolutionless) mas- 
ter equation ([9]) since besides the Born approximation, 
the (first) Markovian approximation of replacing p(t') 
with p(t) in Eq. (UJ) seems to be an additional approx- 
imation made on the time-local master equation. But it 
has been show n 43 ' 46 ' 48 ' 51 that this may not be the case. 
In many example s 43 ' 46 ' 48 ' 51 , the time-convolutionless ap- 
proach works better than the time-nonlocal approach 
when the exact dynamics is available and is used to test 
the perturbative non-Markovian theory based on these 
two approaches. The Markovian approximation that we 
refer to here corresponds to the (second) Markovian ap- 
proximation where the bath (reservoir) correlation func- 
tions are <5-correlatcd in time. In this (second) Markovian 
limit, one may change in Eq. © the integration variable 
t' to t — t— t' and then extend the upper limit of the time 
t integral to infinity (i.e., t — > oo) as the bath correla- 
tion functions (kernels) are ^-correlated in time and thus 
sharply peak at the lower limit r = of the integral. Re- 
cently, there are several investigations^— of construct- 
ing different measures of non-Markovianity to quantify 
the degree of non-Markovian behavior of the quantum 
time evolutions of general systems in contact with an en- 
vironment. In this article, we do not concern ourselves 
with determining the degree of non-Markovian charac- 
ter as investigated in Refs. Io7l - l70l The non-Markovian 
process here refers to that we do not make the (second) 
Markovian approximation of assuming the bath (reser- 
voir) correlation functions being <5-correlated in time on 
Eq. © to obtain the second-order time-convolutionless 
master equation so the resultant quantum dynamics of 
the system state is not Markovian. In other words, the 
influence of the coarse-grained environment causes non- 
local noise correlations and the memory effects of the 
non-Markovian dynamical process are contained in the 
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time- dependent decoherence, damping and/or diffusion 
coefficients of the time-convolutionless master equation 
rather than time-independent ones in the Markovian case. 

We will derive conditional number-resolved (or n- 
resolved) and unconditional quantum master equations 
for the reduced density matrix of the NMR up to sec- 
ond order in the effective tunneling amplitude and in the 
NMR-thermal-bath coupling strength. The n-resolved 
master equatio n 20 i 21 i 24 ~ 27 i 29 i 31 i 32 describes the dynamics 
of the reduced NMR system state conditioned on the 
number n of the electrons that have tunneled through 
the QPC detector in time interval of (0,t), and is thus 
ready to be used to calculate the transport properties, 
such as the transport current. The unconditional master 
equation can be obtained by summing all possible num- 
bers of electrons n in the right lead (drain) of the QPC. 
We will present the derivation of the non-Markovian 
(time-convolutionless form) master equation of the re- 
duced density matrix of the NMR system by consider- 
ing only the non-equilibrium QPC fermionic reservoirs 
first, and will include the effect of the equilibrium ther- 
mal bosonic bath into the derived master equation later. 
To proceed the derivation, it is convenient to go to the 
interaction pictur o 48 ' 65 with respect to Hq = Hs+H\ ca( i s . 
The dynamics of the entire system is determined by the 
time-dependent tunneling Hamiltonian in the interaction 
picture 

Hi(t) = -fftunneling(i) 

= E ^ + Wit)} e< € °- e ?> /h cl tk c D>q + ^.iq) 

ft . q 

where x(t) — x cos(uj t) 4- (p/muj ) sin(u; £). By rewriting 
x(t) 



2 2muj 



2 2muj 



the interaction (tunneling) Hamiltonian, Eq. (jlOl) . can 
be written in the form of 

Hi(t) = E S kq (t)Ft q (t) + Sl q (t)F kq (t), (12) 

k,q 

where 

F kq {t) = e-^-O^cs^ (13) 

is the reservoir operator and 

S kq (t) = [Pi + e iuot P 2 + e'^ ot P 3 ] (14) 

is an operator in a discrete Fourier decomposition 30 act- 
ing on the Hilbert space of the NMR system. Intr oducing 
the dimensionless characteristic length xq = \Jhjmuj 
and momentum po = \/mhuj , we may write 

Pi = T kq , 

P2 = Xkq ( 



X 


i P 


2x 




X 




2x 





(15) 
(16) 

(17) 



where Xkq = Xkq%o has a dimension the same as T kq . The 
form of Eq. ([14]) indicates that there are three different 
jump processes that involve no excitation, excitation, and 
relaxation of the energy quanta of the NMR, respectively. 
Pi is associated with elastic (no excitation) tunneling of 
electrons through the QPC junction and P2 (P3) is asso- 
ciated with inelastic excitation (relaxation) of electrons 
tunneling through the QPC with an energy transfer fwj . 
The energy which relaxes (excites) in response is provided 
by the NMR. By regarding the tunneling Hamiltonian, 
Eq. p2p . as a perturbative interaction Hamiltonian, one 
can see that the first-order term vanishes, i.e., Eq. © 
is satisfied, for the density matrix of the QPC reservoirs 
(leads) given by R = pf cad ® p£ ad , where 



(H. 



Plead 



Tr^e-^ead-^)] 



Here iVj = J2 k c l k Cl , k > Ms and f-D are the chemical po- 
tentials which determine the applied QPC bias voltage, 
eV = [is — /AD, and (3 = 1/ (fc^T) is the inverse temper- 
ature. One may then obtain the second-order (Born ap- 
proximation) time-convolutionless non-Markovian mas- 
ter equation for the reduced density matrix of the NMR 
system by substituting Eq. (fl"2j) into Eq. ([9]). However, 
in order to make contact of the NMR system with the 
QPC detector output current, it will be convenient to 
obtain an n- resolved master equatio n 20 ' 21 i 24 ~ 27 i 29 ' 31 i 32 . 
The rotating wave approximation and Markovian approx- 
imation are known to be pretty good approximations 
for systems in quantum optics. The rotating wave ap- 
proximation is a good approximation provided that the 
strength of the dissipative corrections or the relaxation 
rate, denoted generically as 7^, is small compared to the 
minimum nonzero system frequency difference (energy 
difference/ft) involved in the problem. In the present 
case, this implies that the generic rate 7^ in conditional 
n-resolved and unconditional master equations should 
satisfy the condition of 77? -C wq. The relevant phys- 
ical condition for the Markovian approximation is that 
the bath correlation time is very small compared to the 
typical system response time. But since the (renormal- 
ized) resonant frequency u>q of a NMR is typically in 
the range of a few hundred KHz to a few GHz which 
is much smaller than typical optical frequency of 10 15 Hz 
and since the reservoir correlation time in solid-state de- 
vices may not be much shorter than the typical system 
response time, we will not make the Markovian approx- 
imation and the pre-trace and post-trace rotating wave 



l = S,D. 



(18) 



approximation: 



.71-73 



in our derivation of the master equa- 



tion. By first identifying the jump operator terms and 
partially taking trace over the microscopic degrees of free- 
dom of the QPC reservoirs but keeping track of the num- 
ber n of electrons that have tunneled through the QPC 
detector during the time period (0,i), and then chang- 
ing from the interaction picture to the Schrodinger pic- 
ture, we can obtain from Eqs. ([§]) and <\12\ the time- 
convolutionless non-Markovian n-resolved (conditional) 
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master equation as 



i 



i 



dU 



E 



x 



k,q;k' ,q' 



| f ; t ',(Mi)[4,w^,,(ti)pr w 

-Sj v (*i)p£ +1) (t)Sfc, (*) 

+*4 n) (t)Sl q ,(h)S kq (t)'_ 
+F k * q . k , q ,(t,t 1 )[s kq (t)S^(tOp^ (t) 
-5t, g ,(ii)pk" +1) (*)5fe,(t) 



+£damp [P/j^ (*)] 



(19) 



where we have also included the intrinsic dissipation ef- 
fect of the NMR induced by interacting with a non- 
Markovian thermal bosonic environment in the last 
term of Eq. (|19j) . The mode-dependent symmetric and 
anti-symmetric two-time reservoir correlation function, 

F fe, g; fc',g'(*>*i) and ^fc,g;fe',8'(*>*i) in E q- 0S> are respec- 
tively 



Ft 



kq;k' q' 



(1 - N Dq ) + (1 - 7V Sfe ) jv d J 



and 

Fk,q;k',q' (Ml) 



xe 

1 

: 2 
1 



i(e*-<£)(*-ti)/*<$ 



k,q;k' ,q' ) 



(20) 



{JV Sfc (l-JVD,)-(l-JV Sfc )JVD,} 



xei^-eHMO/^^.^^. (21) 

Here the notation (...) indicates the expectation 
value over the initial density matrix of the reser- 
voirs and consequently F^ q . k , q ,(t,ti) and F£ q . k , ql (t,ti) 
are given by the combination of the Fermi distribu- 



and Nu q 



t inn functions 7V S fe = e /3 ( e ^" /is ) + 1 
e ^( e g -md) _j_ 2 f ^ ne i e ft (source) and right (drain) 

reservoirs of the QPC^l 3 ^. 

By taking into account relevant tunneling amplitudes 
and summing over the wave vectors of the QPC reser- 
voirs, the structure of the influence of the QPC reservoirs 
on the dynamics of the NMR system is then characterized 
by the symmetric and anti-symmetric two-time reser- 
voir correlation kernels, £)*,,,*' ,,' A l, q B k, q Fk, mk ' , q > (*i *i) 



and T,k, q ,k', q ' A l, q B k,q F k, q ;k>, q ' (Ml), where the value 
of A k<q and B k . q could be any one of the tunneling 

amplitudes, T k>q and Xk,q — Xk,q%o- In the con- 
tinuous limit, the summation of the QPC reservoir 
modes can be replaced by the continuous integrations, 
EfcE 9 -> I S d4<kq 9L{4)9R{tq ), wh ere the energy- 
dependent densities of states gs(4) and 9D{t q ) are 
introduced for left and right QPC electron reservoirs, 
respectively. In principle, the tunneling amplitudes, 
T k ,q = T (ef , ) and x M = x (ef , ) , are also energy- 
dependent. We may deal with any realistic energy func- 
tion form of the densities of states and tunneling am- 
plitudes to take into account the memory effect of the 
QPC reservoir on the electron transport and the NMR 
system in our non-Markovian treatment. For simplic- 
ity, we follow several non-Markovian electron transport 
studie a 33 i 36 ~ 39 ' 41 ~ 43 ' 74 ' 75 by considering a spectral density 
with energy-dependent densities of states and tunneling 
amplitudes absorbed in a Lorentzian form as 

J A , B {el ef ) = ^ (4, ef) B (4, e?) 9l(4 )te(e?) 



A< B 00 glg° R A 2 



A 



(22) 



where the cut-off energy A e characterizes the width of the 
Lorentzian energy-dependent distribution, the parameter 
Ei denotes the effect of the variation of the QPC junc- 
tion barrier potential^ 3 , due to the interaction with the 
NMR, Aqo, -Booj g\ and g R are energy-independent tun- 
neling amplitudes and densities of states near the aver- 
age chemical potential. Physically, this spectral density 
of Eq. (f!?2j) means that given an electron state with a 
fixed energy e k in the left lead, the electron can tunnel 
into the electron energy states of the right lead with a 
central energy at + Ei and a Lorentzian width A e . 
In the limit of A e — > and in the absence of the inter- 
action with the NMR (i.e., Ei = 0), the QPC spectral 
density, Eq. (1221) . is proportional to 8 (ef — e q ) that rep- 
resents the resonant tunneling process. In the opposite 
case of the cut-off energy A e — > oo, the QPC spectral den- 
sity, Eq. (|2"2"|) . becomes energy- independent and reduces 
to the constant WBL spectral density used in the liter- 
ature. The average (effective) zero-temperature electron 
tunneling conductances (G/e 2 ) through the QPC bar- 
rier in the WBL can be written as (27r/ft)Aj £?ooSL3 l R- 
Compared with the energy-dependent spectral density, 
the WBL one that assumes all electron states in the left 
reservoir having equal likelihood to tunnel to all the elec- 
tron states in the right reservoir regardless their energies 
may not be a very good physical approximation after all. 

We note that the dynamical behaviors of the NMR- 
QPC system are sensitive to the actual energy depen- 
dence and the bandwidth of the QPC spectral density 
(which may be different from the simple Lorentzian form 
considered here). The realistic energy dependence or 
function form of the densities of the states and the tun- 
neling amplitudes in the spectral density depends on the 
detailed QPC electronic structure. Here we perform a 
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model calculation for the QPC-NMR system using a sim- 
ple Lorentzian spectral density to study the influence of 
finite bandwidth (cut-off energy) and memory effect on 
the NMR system dynamics. We will show later in our 
numerical treatment of the non-Markovian NMR-QPC 
system that for the parameters we choose, when the 
bandwidth of the Lorentzian spectral density of Eq. (|22|) 
is about A e < bhuj 0l the time-dependent coefficients, 
the dynamical variables of the NMR and the currents 
through the QPC differ significantly from their Marko- 
vian WBL counterparts. This can be understood as fol- 
lows. As discussed earlier, in the limit of A e — > 0, only 
one channel ef — = Ei is involved in the electron tun- 
neling processes across the QPC barrier. The opposite 
limit of a very large bandwidth, A e » | ef — — E% | , then 
leads to a channel-mixture regime^ where great portions 
of all possible ef ^ transitions of electron tunnel- 
ing between the source (left reservoir or lead) and the 
drain (right reservoir or lead) are allowed with weight de- 
termined by Ja,b ( e k' e ^) an d with randomness coming 
from electron scattering determined by the band struc- 
ture associated with the geometry of the metallic gates 
in the QPC. The electron-tunneling processes with more 
random channel mixture will reduce the QPC reservoir 
correlation time^ and therefore suppress the QPC reser- 
voir memory effect on the NMR dynamics. Thus the 
non-Markovian processes will become significant if the 
QPC electronic structure can be designed or engineered 
to have a spectral density with a narrow bandwidth com- 
parable to the (renormalized) resonant frequency of the 
NMR as the QPC reservoir correlation time in this case 
is comparable to the NMR system response time (see also 
the discussions regarding Figs. [2] and [3] in Sec. IVII|) . The 
typical frequency of NMR is in the range of a few hundred 
KHz to a few GHz. Thus the condition of being able to 
observe a significant non-Markovian finite-bandwidth be- 
havior of A e < 5hajo suggests that the bandwidth of the 
QPC spectral density should be in the range of about 
1 - 20 fieV. 

With the specified spectral density Eq. (|2"2"j) and the 
help of Eqs. (|20|) and (f2Tj) . one can rewrite the n-resolved 
master equation (|19p into the following form 



P% ) (t) = -j_[H s , P %>(t)} 

+ 9 ^{j^ft{t,eV+^ i ) 

3 
i=l 

+ £damp[,4 n) (i)]. (23) 

Here P l is defined in Eqs. (HU) dTTJ) , P = ELi P * = 



Pi + P2 + P3, the values of the frequency LOi are given by 
uji = 0, UJ2 = — UJ3 = uj 0i and H.c. denotes the Hermitian 
conjugate of all the previous terms in the curly bracket of 
Eq. (|2"3")l . By the change of the new variables, ujf = ef — 
Us and u>® = — /xd, the time-dependent coefficients 
/f(s) m P^l) can be written as following forms: 



f+(t,eV)= [f F (t,eV)Y 



= - dr 



00 /»oo 



OO J — OO 



A 2 



1 



e^t + 1 



f - 



xe 



+i(u,S-u;°+eV)T/h 



(24) 



f+(t,-eV)= [f B (t,-eV)] 1 



OO />OG 



— OO «/ —OO 



A 2 



(ul-ojV + eV-E^+Al 



x 1 



f 



1 



e^l + lj e P"? +1 

xe -H(«f-o,f+eV)r/ft j (25) 



/+ (i, eV ± hw a ) = [f F (t, eV ± hu )] f 



S Jo 



dr 



OO />00 



duldujf 



— OO J — OO 



A 2 



{LoS-uD + eV-Etf+kl 



e P< + l 



eP»? + 1 



xe +t(w£-w£+eV±fiui )T/ft ^6) 



and 



/+ (t, -eV T nw ) = [f B (t, -eV T * 



OO /*00 



OO J — OO 



A 2 



(ul-^+eV-E^+Al 



1 



x- 



e^< + 1 

Physically, f F (t, eV) and f B (t, —eV) describe the mem- 
ory effects on the NMR system induced by the elastic 
electron tunneling processes in the QPC reservoirs with 
no excitation of the NMR. The time-dependent coeffi- 
cients fp (t, eV ± huj ) and f B (t, —eV =p fc^o) describe 



f - 



1 



e m + 1 
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the memory effects on the NMR system caused by the in- 
elastic electron tunneling processes that involve the NMR 
excitation and relaxation, respectively. 

The effect of the thermal bosonic environment on the 
master equation in the last line of Eq. (f2l?f can be derived 
also up to second order in system-environment coupling 
strength^ 8 - - — and the result is given as 



£dampbi*°(*)] = - r 

i 

1 

~W 
1 



7o (t) [x,{ P , P ^(t)} 

D (t) 



■r. 



h (t) 



where the timc-dcpcndcnt coefficients are 



2 
m 



dr cos(w T)ry(r), 



(28) 



(29) 



where A D is the cut-off frequency, 7 is a constant char- 
acterizing the strength of the interaction with the envi- 
ronment, and m is the mass of the NMR. For simplicity, 
we will take the commonly used spectral density of an 
Ohmic bath, i.e., n = 1 case in Eq. (f3"B")) . 

The ra-resolved master equation (|23|) with Eq. (|28|) was 
derived without making the Markovian and the pre-trace 
and post-trace rotating wave approximations, and the 
only approximations we use are the second-order pertur- 
bation theory, the initially factorized system-bath den- 
sity matrix, and the forms of the spectral densities of 
Eqs. (|2"2"|) and (|36[) . So the n-resolved master equation 
is valid for arbitrary bias voltages and environment of 
temperature, as long as the perturbation theory that we 
use up to second order in the system-QPC and system- 
environment coupling strength holds. 



IV. UNCONDITIONAL MASTER EQUATION 
AND MARKOVIAN LIMIT 



7o(t) 



moJo Jq 



dr sin(w r)?7(r), 



D (t) = h / dr cos{u t)v{t), 
Jo 



h (t) = -- 



muj J 



dr sin(w r)j/(r). 



(30) 



(31) 



(32) 



Here Cl^(t) is the frequency shift due to the coupling to 
the thermal environment, 7o (t) is the dissipative coef- 
ficient, and D {t) and h (t) represent the diffusion co- 
efficients. The two kernels t)(t) and v{t) appearing in 
Eqs. (f29"|) — (|32|) are so-call dissipation and noise kernels , 
respectively, and are defined as 



t?(t) 



duiJ (u) sin (lot) , 



(33) 



and 



divJ (lj) cos (lot) coth(/3hw/ 2) , (34) 



where J{u>) is the spectral density of the bosonic envi- 
ronment defined as 



J(u) 



^ 2m„w„ 



Again, we could, in principle, deal with any given form 
of the spectral density. But as a particular example, we 
use the following form of spectral density with a Lorentz- 
Drude cut-off function to specify the environment ^ 50 ' 58 



A. Unconditional master equation 

In this subsection, we present the unconditional master 
equation for the reduced density matrix of the NMR sys- 
tem. Statistically, the unconditional master equation can 
be straightforwardly obtained by summing up Eq. (1231) 

over all possible electron number n, i.e., p R (t) = J2 n Pr '■ 
Despite the different nature between the non-equilibrium 
fermionic QPC reservoir and the thermal bosonic en- 
vironment, by combining the relevant terms together, 
the unconditional non-Markovian master equation can 
be cast into a simple form similar to the non-Markovian 
quantum Brownian motion master equation as 



PR(t) 



H. 



sys 



^(ul{t) + & {t))x\p R {t) 



-j:he(t)+lo(t))[x,{p,p R (t)}] 
-±(D e (t)+D (t))[x, [x,p R (t)]] 
+ jp(he(t) + h (t))[x,\p,p R {t)}}. 



(37) 



The whole non-Markovian character of the dynamics of 
the NMR system is contained in the time-dependent co- 
efficients appearing in the master equation. The time- 
dependent coefficients that come from the QPC electron 
reservoirs are denoted with a subscript e. The frequency 
(35) renormalization ui^(i), the damping coefficient 7 e (£), the 
decoherence coefficient D e (t) and the diffusion coefficient 
h e (t) are, respectively, given by 



(38) 



J(a,) = -m 7 u,(£) 



At 



(36) 



7«(*) 



2itmuJo 



Re Kf (*)-&*(*)]» 



(39) 
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D e (t) 



h e (t) 



h 2 G x 
2tt 



L Re [£(*) + &'(*)], (40) 



h 2 r 

2-Kmui 



Im^W-^W], (41) 



where 





2tt 


s£a& 1. 


1 2 






(42) 






(t,eV 


+ ftwo) + 


ft 


(t,-eV- 


foj ), (43) 


£?(*) 




(i,eV 


+ fejo) - 




(f,-eV- 


fiwo), (44) 






(i,eV 


— huj ) + 


f + B 


(t,-eV + 


fiuo), (45) 


% (*) 




(t,eV 


- hw ) - 


n 


(t,-eV + 


Bw B ). (46) 



upper limit t of the time r integrals in Eqs. (|24 p -(|27 p 
for QPC reservoirs and in Eqs. (I29|) - (l32t for the thermal 
bosonic bath can be taken to t — > oo. Another com- 
monly used assumption in the Markovian limit is the so- 
called WBL approximation. This assumption may not 
seem essential to evaluate the integrations if one already 
makes the very short bath correlation time approxima- 
tion. But the assumption of very short bath correlation 
times can be justified in various models of the bath spec- 
tral densities with very large cut-off energies. This has 
been demonstrated for the bosonic Ohmic bathe r 50 1 58 . 
We will show in Sec. I VIII that this is also the case for 
the simple Lorentzian spectral density of Eq. (1221) for the 
non-equilibrium QPC fermionic reservoirs. Thus, if we 
take the Markovian approximation of very short corre- 
lation times (integration limit t — > oo), then the time- 
dependent coefficients in the master equations (|23|) . (l28l) 
and (|37[) become time-independent. Specifically, using 
the relation 



We note that due to the interaction with the ther- 
mal environment, the frequency shift term Ctg(t) in Eq. 
(|37| diverges as the cut-off frequency A G — > oo and 
thus is not physical^ 8 - -50 ' 58 . Therefore, a regularization 
procedur o 76 ' 77 is needed for the frequency renormaliza- „t 

tion. We adopt the view of the renormalizatio n 76 i 77 to re- lim / dre 1 ^^""^ = ir5 (oj — ujq) + iPV 
gard the frequency in the original Hamiltonian as a finite t^oo j Q 
renormalized frequency u> and add a frequency counter 
term with a frequency^ 8 - -50 ' 58 defined as 



1 



(48) 







1 

rn 



E 



A 



m n w 2 



2 

m jo 



duj 



J(w) 



(47) 



where PV indicates the Cauchy principle value, the co- 
efficients /± (B) in Eqs. (IMl)- ([27)l coming from the QPC 
reservoirs can be written as 

lim (i, y) = W F[B) [y] ± iQ F{B) [y] , (49) 



to cancel at large times the frequency shift ^(t). Similar 
to this reasoning, another frequency counter term with 
a frequency io\ is introduced 21 to compensate at long 
times the frequency shift w 2 (i) induced as a result of 
the coupling to the QPC leads (reservoirs). The physical 
frequency in this case is then w 2 (i) = w 2 + o> 2 (i) + f2 2 (t) + 
u; 2 + f2 2 and approaches the finite renormalized frequency 
lo at large times. The non-Markovian master equations 
(f23"l) and ([37jl are the main results of this paper. 



where 



W F [y] = 7T 
W B [y] = 7T 
e F [y] -- 



A 2 



1 _ e -y/k B T ( e V-y + E t ) 2 + A 2 



y 



A 2 



B. Markovian limit 

Next, we show that by taking appropriate limits, our 
non-Markovian master equations can recover the various 
Markovian master equations reported in the literature. 
The Markovian approximation is valid when the bath 
correlation time is much smaller than the characteristic 
time scale of the system of interest. The bath correla- 
tion time is determined by the bath correlation kernels 
(functions) and depends on the form of the bath spectral 
density. We will perform the numerical calculation of the 
QPC reservoir correlation time in the next section to in- 
vestigate how the reservoir correlation time is varied as 
a function of various parameters in the problem. 

Here, if we nevertheless take the Markovian approxi- 
mation of very short bath correlation times of the QPC 
reservoirs and of the thermal bosonic environment, this 
is equivalent to assuming that the bath correlation func- 
tions (kernels) are S— correlated in time and thus the 



1 _ e -v/kBT ( e y + y + E .f + A 2 

A 2 



, (50) 



(51) 



OO fOC 



oo «/ — oo 

1 



dw^duj® — s ; 

1 



1 - 



( 1 

xPV 



- + eV + y r 



(52) 



s dio D 



A 2 



/OO nOO 




- uf + eV - y 



(53) 



We note here that we separate the Markovian ap- 
proximation from the WBL approximation although 
the Markovian approximation can often be justified 
by considering a very large cut-off energy. So the 
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cut-off energy A e remains in Eqs. (|50|) -(f55 )) . If the 
WBL (A e — > oo) is taken, then the real parts 
of fp(B) ~ ^ 00 ' ?/) wnen multiplied by the factor 
29'x,3flj4oo-Boo/fi become the finite-temperature forward 
(backward) Markovian WBL electron tunneling rates of 
rf£(y) = (2n/h)g L g R A oB o/[l-exp(-y/k B T)} in the 
literature 1 ^ 21 i 25 i 28 where the value of A 00 and B o 
could be either one of the tunneling amplitudes Too and 
Xoo- As a result, in the Markovian limit, the frequency 
renormalization, the damping coefficient, the decoher- 
ence coefficient, and the diffusion coefficient in Eqs. 
(|4T1) due to the QPC reservoirs in the WBL (A e - 
become, respectively 



00) 



hG x 



-{Q F [eV + hw ] -S B [~eV- 
+e F [eV-fiu ] -e B [-eVA 



M 



D 



M 



h -G 

m 

II \JTxx 



(55) 



(eV + huj ) coth 
+ (eV — hu! ) coth 



eV + Hlu 



2k B T 
eV — hu 



hG 3 



2k B T 

■ {Q F [eV + hoj ] +e B [-eV - hu>o] 



(56) 



2ttitiuj 

-&f [eV - hw ] -Q B [-eV + huj ]\57) 

Similarly, the frequency renormalization, the damping 
coefficient and the diffusion coefficients in Eqs. (|2"9" ]) -([3lZ ]) 
due to the Ohmic thermal environment in the Markovian 
WBL (A Q — > 00) can also be obtained as 48 



A 



= / duj 



J(u) 



7o M 



D 



M 



7, 



myhug coth 



huj 

2k B T 



h? = jk B T 



n=— o 



(58) 
(59) 

(60) 
(61) 



where v n = 2imk B T '/H are known as the Matsubara fre- 
quencies. 

We note that our unconditional non-Markovian mas- 
ter equation (f3"7| in the Markovian WBL recovers the 
Markovian master equation (2.12) in Ref. |2jJ. In the 
special case of the zero-temperature and high-bias limit 
where D^f = mj^' 1 eV, we recover Eq. (6) of Ref. [lj^ 
if the coefficients coming from the contributions of the 
thermal bosonic bath are neglected. Similarly, the con- 
ditional non-Markovian n-resolved master equation (|23[) 
also reduces to the Markovian n-resolved master equation 
(2.9) in Ref. Hi]. Taking the high-temperature limit of the 
bosonic environment and Fourier-transforming in the n 
index, the conditional n-resolved master equation (|23[) in 
the Markovian limit also reduces to Eq. (2) of Ref. 



if the dc bias case is considered and the transmission 
phase 77 is set to zero in Ref. 20, Again, considering only 
the QPC reservoirs and in the special case of the zero- 
temperature and high-bias limits, one obtains Eq. (5) of 
Ref. [HI in the Markovian limit from the conditional n- 
resolved master equation (f23|). 



V. DYNAMICS OF THE NMR 



Using the master equation (|37j) . we can obtain the 
equation of motion for the mean or expectation value 
of any physical operation O of the NMR by calculating 
= Tr[Opn(t)]. So the equations of motion of the 
mean (expectation value) of the position and the mo- 
mentum are 



d(x(t)) _ (p(t)) 

dt m 
d(p(t)) 



dt 



-muj 2 p (t) (x(t)) 



-2[7e(*) + 7o(*)]<P(*)>» 

and for the second moments we obtain 
d(x 2 (t)) 1 



d{p 2 («)) 



= -mu 2 Jt) ({x, p}(t)) 



dt "" py 

-4[ 7e (i)+ 7o (i)] (P 2 (*)) 
+2[D e (t)+D (t)], 
d({x,p}(t)) _(p 2 (t)) 



(62) 



(63) 



(64) 



(65) 



dt 



= 2- 



2m^ p {t) (x 2 (t)) 



-2[ 7e (i)+ 7o (i)] ({x, p}(t)) 

+2 [h e (t) + h a (t)} . (66) 



Combining Eqs. ([62]) and (|63|) yields 

d(x(t)) 



d 2 (x(t)) 
dt 2 



+ 2 7toi (t) 



dt 



«4(t)(x(t)) = 0, (67) 



where j to t = 7e + 7o- 

One may in principle solve the time evolutions of the 
differential equations (|62"T) -(|66 [) and the numerical results 
will be presented in Sec. IVIII Simple analytical expres- 
sions of the steady-state (t — s- 00) solutions can, however, 
be obtained asZ 8 - 



(P 2 ) t 



(P)t- 
0, 



0. 



(68) 
(69) 



lim 



lim 



D e (t) + D (t) 
m 2mui 2 (t) \m[y e (t) + 7 G (t)] 

+2 [h e (t) + h {t)\ 
D e {t)+D {t) 



*-"»2[ 7e (t)+ 7o (t)] 



(70) 
(71) 
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For the moment, let us consider the case where the in- 
fluence of the thermal environment is neglected. We also 
note that, for typical values of finite electric reservoir 
temperatures and finite electric bias voltages, the diffu- 
sion coefficient h e (t)/h is generally much smaller than 
D e (t)/po and J e {t) and thus is often neglected. In this 
case, we obtain from Eqs. ([70]) . (joo) and ([So) the steady- 
state (a; 2 )t^oo in the WBL as 



1 



(eV + FiLJo) coth 



eV + hio 



+ (eV — Hoj ) coth 



2k B T 
eV — hw 



2k B T 



(72) 



At zero temperature (k B T 
hjj n ), we have (a: 2 ) ~- 



0) and low voltages (eV -C 
. In this case, the NMR is 
in the ground state and is independent of the bias volt- 
age as the bias voltage is unable to excite the NMR from 
its ground state. On the other hand, at high voltages 
(eV hu ), the NMR is no longer in the ground state 
and (a; 2 ) w At high temperatures (k B T > eV, 

huj ), the quantum mean square of the position of the 
NMR becomes (x 2 ) s» ss^j. that is expected from a clas- 
sical oscillator in thermal equilibrium. 



VI. TRANSPORT CURRENT 

With the n-resolved time-convolutionless master equa- 
tion for p^\t), one is readily able to compute transport 
current I (t) = e^^-, where (N (t)) = £ n nP (n. *) = 
^2 n nTr p R (i) is the expectation value of the num- 
ber of electrons that have tunneled into the right lead 
(drain) in time t. Here Tr means tracing the density 
matrix over the degrees of freedom of the NMR system. 

Inserting Eq. (j23|) into I (t) = e J2 n "Tr p R (t) gives 
rise t^2^-^L2M^l 



I it) 9l9°n 



£ Tr { [/+ (t, eV + hui) P}Ppn{t) 



-fZ(t,-eV-nui)PP}p R (t) 



H.c 



}.(73) 



Here p R (t) — ^2 n P%\t) i s tne unconditional density 
matrix of the NMR. This non-Markovian average cur- 
rent is valid for arbitrary QPC lead temperatures and 
arbitrary bias voltages as long as the second-order per- 
turbation theory holds. We follow Ref. Hi] to categorize 
the non-Markovian average current, Eq. (|73p. into four 
physically distinct contributions. Using the definition of 
p = E<=i p i and thc definition of P l in Eqs. ([ ^ -(fl? ) . 
we write the non-Markovian average current, Eq. (|73p . 
as 



J(t) _ 

-^position 

e e 



(t) , I P (t) I QM (t) I{ X ,P}(t) 



(74) 



The first term in Eq. (f74"|) depends on the oscillation po- 
sition of the NMR and can be written as 



'position 



(*) 



1 

2^ 



Ai (t) [G + G x (x(t))] 



47T 



[A 2 (t) + A 3 (*)] 



x [G.^^+G^^fi))] . (75) 

It reduces in the Markovian WBL to the so-called Ohmic- 
like part of the current proportional to the conductance 
defined in Ref. 21. In Eq. (|75|) . thc term with conduc- 
tance Go = ^-g L g R \T 00 \ 2 represents the current through 
the isolated QPC junction, and the remaining terms due 
to the coupling to the NMR with conductances G x = 
%&g R MT$oXoo] and G xx = f glg R \ X m\ 2 contribute 
to the nonlinear part of current- voltage characteristi c 21 ' 79 
as the state of the NMR will depend on the bias volt- 
age. The time-dependent coefficients Aj(t) in the non- 
Markovian region can be written as 

A 1 (t) = f+(t,eV)-f+(t,-eV) 
+f F (t,eV)~f B (t,-eV) 
= 2Rc [/+ (t, eV) - /+ (t, -eV)) , (76) 
A 2 (t) = /+ (t, eV + hio ) - /+ (t, -eV - hio ) 
+f F (t, eV + Hlj ) - f B (t, -eV - hw ) 
= 2Re [£(*)], (77) 
A 3 (f) = # (t, - ftw ) - /+ (t, -eV + 

+/ F (t, el/ - fiw„) - f B (t, -eV + Tiwo) 
= 2Re[£ 2 a (t)], (78) 

where £f (i) and £f (t) are defined in Eqs. (jSJl and (|46|) . 
respectively. The second term in Eq. (|74")l . Ip(t)/e, is 
proportional to the average velocity of the oscillator 



J P (t) 



G 



^_[A 2 (t)-A 3 (t)]^(f)). (79) 



This term is nonvanishing only for an asymmetric 
junction 21 , G p = ^g L g R lm[T^ X oo} ^ 0. The third 
term in Eq. (|74p . Jqm (t), derived from the commuta- 
tor of position and momentum operators, is referred to 
as the quantum correction to the currents! 



(*) 



KG, 



8irmuj 



[A 4 (t) - A 5 (t)} . 



(80) 



Here the time-dependent coefficients A 4 (t) and As(i) can 
be written as 

A 4 (t) = (t, + fta; ) + /+ (t, -eV - fiwo) 
+/ F (t, eV + fiw ) + /g (t, -eV - 
= 2Re[ff(t)], (81) 
A 5 (i) = ft («, - fiw ) + /b (*, -eV + b„) 
+/ F (t, el/ - tiuo) + f B (t, -eV + fiw ) 
= 2Re[e|(t)], (82) 
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where £f(t) and £J(i) are denned in Eqs. (g3J| and ([15]). 
respectively. By using Eqs. (l49l) - (|5Tj) in the Markovian 
WBL (i.e., taking A e — > oo), the above three terms that 
contribute to the total current recover their respective 
Markovian versions presented in Ref. I2D : 



rM 

position 



tM 
1 P 



tM 
QM 



eV [Go + 2G x (x(t) + G xx (x 2 (t))] , (83) 



^<p(i)) 

TO 

hG. r 



4mw„ 



(ey + fouj ) coth 



(eV — hid ) coth 



eU + 7kj 
2k B T 
eV — fiw 



2fc B T 



(84) 



(85) 



The last term in Eq. (|74l) . T ry P y(t)/e, is a new term 
that was not discussed in Ref. l2lL This term originates 
from symmetrized product of the position and momen- 
tum operators and can be written as 



L {X,P 



}(t) 



iG x 



4to,Cl! 7T 



[A 6 (t)- A 7 (t)]((xp + px) (t)> 



where 



(86) 



A 6 (t) = ft (t, eV + frw ) - f+ (t, -eV - hw ) 
fp (t, eV + frw ) + f B (t, -eV - hw a ) 
= 2*Tm[ff(t)], (87) 
A 7 (t) = /+ (t, eF - fi^ ) - /+ (*, -eV + fi^o) 
-fp (t, eV - hw ) + f~ (t, -eV + hw ) 
= 2ilm[C 2 {t)], (88) 

and £f(t) and ££(f) are defined in Eqs. flU) and gB]), 
respectively. In the Markovian limit, this term becomes 



{X,P} _ (jxx 



2mU! n TT 



F (eV + hw ) - B (-eV - Hlu ) 

-Of (eV - huj ) + B (~eV + huj 
x((xp + px)(t)), (89) 

where the functions 6^(0;) and Ob (x) are defined in 
Eqs. (|52t and (|53)l . respectively. This extra current term, 
Eq. (|86p. with a coefficient coming from the combina- 
tion of the imaginary parts of the QPC reservoir cor- 
relation functions, was generally ignored in the studies 
of the same problem in the literature^—. The contri- 
bution to the diffusion coefficient h e (t), Eq. (j4"Tj). also 
comes from the imaginary parts of the QPC reservoir 
correlation functions but with a different combination. 
Unlike the diffusion coefficient h e {t) which is generally 
much smaller than the other decoherence and damping 
coefficients for typical parameters and is thus often ne- 
glected, we will show in the next section that Eq. (|86|) 
has a substantial contribution to the total transient cur- 
rent in the non-Markovian case and differs qualitatively 
and quantitatively from its Markovian WBL counterpart, 



(a) 
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% 

* 5 




0.05 0.1 



v.. 



(b) 



eV 
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eV 
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FIG. 2. (Color online) Real part of (a) the symmetric reser- 
voir correlation kernel (function) F s (t) and (b) the antisym- 
metric reservoir correlation kernel (function) F a (r) at a small 
Lorentzian cut-off energy of A e = hw for different values of 
the bias voltage: (eV/hw ) = 0.1 (blue dot-dashed lines), and 
3 (red dashed lines). The time is in units of u^ 1 . The insets 
in (a) and (b) are for the cases of (A e /Hui ) = 10 and 100, 
respectively. Other parameter used is ksT — 0.1huj o . 



Eq. (|89[) with A e — > 00. Thus it may serve as a witness 
or signature of finite-bandwidth non-Markovian features 
for the coupled NMR-QPC system. 



VII. NUMERICAL RESULTS AND ANALYSIS 

In our numerical calculations, we first concentrate on 
the case where the influence of the thermal environ- 
ment is neglected. This allows us to address the non- 
Markovian effect coming solely from the QPC reservoirs. 
This case where the effect of the QPC reservoirs dom- 
inates over that of the thermal environment may nev- 
ertheless be justified for a much larger relative coupling 
strength of QPC to the NMR and for typical QPC bias 
voltages and reservoir temperatures. We will discuss in 
subsection IVII Bl the case when the effect of the thermal 
environment is included and is comparable to that of the 
QPC reservoirs. 



A. Effect of only the QPC reservoirs 

Figures and [2Kb) show the dependence of real 

parts of the mode-independent two-time symmetric and 
antisymmetric QPC reservoir correlation kernels (func- 
tions) on the timc-diffcrcncc r = t — t\ for various values 
of the cut-off energy A e . The mode- independent two- 
time symmetric and antisymmetric QPC reservoir corre- 
lation kernels (functions) are evaluated by converting the 
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FIG. 3. (Color online) (a) Real part of the symmetric reservoir 
correlation kernel (function) F s (r) as a function of the time- 
difference t = t — ti and the bias voltage eV. The bias voltage 
is in units of Huj and the time is in units of a;" 1 . Other 
parameters used are A e = huj and ksT = 0.1ftu o . (b) Real 
part of the symmetric reservoir correlation kernel (function) 
F s (t) as a function of the time-difference r = t — ti and 
the QPC lead temperature ksT . The lead temperatures is in 
units of huj and the time is in units of a;" 1 . Other parameters 
used are A e = hw and eV = 0.1ftcj o . 



summations over k, k' and q, q' of the mode- dependent 
correlation functions (kernels) of Eqs. (|2"0|) and ([2~Tj) into 
energy integrations with the Lorentzian spectral density 
given by Eq. (|22p . and then performing the energy inte- 
grations. There are several characteristic time scales in 
this non-Markovian problem. The time scale of the NMR 
is about l/u , the time scale of the energy-dependent 
QPC spectral density is about h/A e , the time scale of 
the applied bias is H/eV, the time scale of the QPC 
reservoir temperature is S/fc^T, the time scale of the 
electron tunneling is about an d the time scales 

of the combinations of the electron tunneling rates are 
Po/^e 1 j l/7e*> ^/^e f - We define the time scale at which 
the profiles of the QPC reservoir two-time correlation ker- 
nels (functions) decay as the QPC bath correlation time, 
tb- We can see from Figs, [^a) and[5Jb) that at a given 
low cut-off energy of A e = Hlo and at a low tempera- 
ture of ksT = 0.lhuj Ol the reservoir two-time correlation 
kernel (function) with a high bias voltage of eV = 3fiw 
(in red dashed lines) has a slightly shorter reservoir cor- 
relation time tb than that with a low bias voltage of 
eV = QAhu! (in blue dot-dashed line). But as indicated 
from the insets of Figs. [2fa) andfSfb), the dependence 
of tb on the bias voltage is much weaker than that on 
the cut-off energy A e . The insets show that the larger 
the cut-off energy A e is, the smaller the bath correlation 
time tb is. Moreover, the bath correlation time in this 
case is about tb ~ h/A e . 

Figure [3Ja) shows the symmetric reservoir correlation 
kernel (function) that depends on bias voltage eV and 
time-difference r = t — t% for a cut-off energy of A e = 
1.0hui o . Although the bath correlation time is affected 
mainly by the value of the cut-off energy A e , one can still 
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FIG. 4. (Color online) Real parts of the time-dependent for- 
ward (backward) inelastic electron tunneling coefficients for 
different values of the QPC bias voltage: (eV/hu) ) = 0.1 (in 
blue solid lines), 3 (in red dot-dashed lines), and 10 (in green 
dashed lines) . The Markovian cases with finite cut-off energies 
are plotted in gray dotted lines. The time is in units of ui~ . 
Other parameters used are ksT — 0.1hu o , A e — 100fta) o . The 
small subplots in (a) and (b) are for the cases of (A e /hu ) — 1 
and 10, respectively. Other parameter used is ksT — 0.1hu o . 



see for the Lorentzian spectral density Eq. (|22j) chosen 
here that when the bias voltage decreases, the bath cor- 
relation time tb increases. Figure [3jb) shows the depen- 
dence of the symmetric reservoir correlation kernel (func- 
tion) on temperature ksT. Similar to that on the bias 
voltage, the correlation time increases when the reservoir 
temperature decreases. The effect of the temperature 
on the bath correlation time seems to be stronger than 
that of the bias voltage. The Markovian approximation 
is valid in the case when the bath correlation time tb is 
much smaller than the typical system response time ts- 
The typical response time of our NMR system is about 
the minimum value of (l/u> , ^-l~^ABiV\l^ei l/7e, h/h e ). 

Figures 21 [5] and [5] show the real parts of some typ- 
ical time-dependent coefficients fp^s) °^ ^ ne ^-resolved 
master equation for different values of the bias voltage 
and temperature, respectively. The cut-off energy, A e , 
is also varied in each subplot. Physically, these coeffi- 
cients, if multiplied by the factor 2g1g%AooBoo / h (where 
the value of Aoq and Boo could be either one of the tun- 
neling amplitudes Too and Xoo), correspond to the time- 
dependent finite-temperature forward (backward) inelas- 
tic QPC electron tunneling rates that accompany with 
the absorption or emission of the NMR energy quanta. 
Compared to their Markovian counterparts that are con- 
stants in time plotted in dashed lines in Figs. |H[5]and[6l 
the memory effects of the QPC reservoirs are contained in 
the time-dependent coefficients. As the cut-off energy A e 
is increased, the profiles (peaks) of the time-dependent 
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FIG. 5. (Color online) Real parts of the time-dependent for- 
ward (backward) inelastic electron tunneling coefficients for 
different values of the QPC bias voltage: (eV/hu) ) = 0.1 (in 
blue solid lines), 3 (in red dot-dashed lines), and 10 (in green 
dashed lines). The Markovian cases with finite cut-off energy 
are plotted in gray dotted lines. The time is in units of w„ ■ 
Other parameters used are fcsT = 0.1ftu o , A e = l00hui o . The 
small subplots in (a) and (b) are for the cases of (A e /hoj ) = 1 
and 10, respectively. Other parameter used is fcsT = 0.1huu o . 

coefficients become higher and the half-maximum widths 
become narrower. Moreover, the positions of the peaks 
of the profiles also shift to the short time region. In 
other words, the memory effects of the non-Markovian 
time-dependent coefficients persist for longer times for 
smaller values of the cut-off energy. We distinguish 
the Markovian case, where the cut-off energy A e is fi- 
nite, from the Markovian WBL case, where the spec- 
tral density becomes energy-independent as A e — > oo. 
So one can also see that the time-dependent coefficients 
approach their respective long-time Markovian counter- 
parts of Eqs. ([50)1 and (ToTT) with finite cut-off energies 
(bandwidths). Furthermore, the non-Markovian coeffi- 
cients with larger cut-off energies approach more rapidly 
in time to their Markovian counterparts. 

Generally speaking, for the inelastic forward tunnel- 
ing coefficients in Figs. Q|a), Ufa), Et a ) and Etb), when 
the values of the temperature and bias voltage are higher, 
the long-time asymptotic (Markovian) values of the time- 
dependent coefficients are larger but the time scales for 
the time-dependent coefficients approaching their long- 
time values become shorter. The inelastic (emission) 
backward tunneling coefficients in Fig. @] (b) approach, 
at a rather low temperature of ksT — 0.1Huj o , approx- 
imately zero (the Markovian value) at long times since 
the argument of (— eV — hw a ) is negative [see Eq. (ToTT) ]. 
Similarly, since the argument of (eV — fiw ) for a low bias 
voltage of eV — 0.1fouj o is negative, the inelastic (emis- 
sion) forward tunneling coefficients in blue solid lines in 
Fig. EK a ) also approach, at a low temperature, approxi- 
mately zero (the Markovian value) at long times. On the 



FIG. 6. (Color online) Real parts of the time-dependent for- 
ward (backward) inelastic electron tunneling coefficients for 
different values of the QPC lead temperatures: {UbT /hu ) = 
0.1 (in blue solid lines), 1 (in red dot-dashed lines), and 10 
(in green dashed lines). The Markovian cases with finite cut- 
off energy are plotted in gray dotted lines. The time is in 
units of u„. Other parameters used are eV = 0.1ftu o , and 
A e = lOOhwo. The inset in each subplot is for the case 
of eV = Q.lhuio and the small Lorentzian cut-off energy of 

A e = hw a . 



other hand in Fig.[5jb), since the argument (— eV + huj ) 
for eV — 0.lhuj o is positive, the phonon-assisted back- 
ward tunneling is allowed even in the negatively biased 
direction at a low temperature and thus the inelastic tun- 
neling coefficients in blue solid lines become finite at large 
times. One can see from Figs. 0] and El that if the bias 
voltage is increased, the inelastic forward tunneling coef- 
ficients in Figs.|U[a) and Eta) are enhanced while the in- 
elastic backward tunneling coefficients in Figs. [4f b) and 
Efb) are suppressed. Unlike the cases of eV = O.lfkv 
in Figs. Bib) and Ufa), the time-dependent inelastic 
tunneling coefficients (in red dot-dashed lines and in 
green dashed lines) in Figs. Etb) and Etc) approach fi- 
nite nonzero values for large temperatures even though 
their arguments of (—eV — ftui ) and (eV — Huj ) are neg- 
ative. Moreover, the non-Markovian coefficients with a 
large cut-off energy of A e = lOOftw shown in Fig. El ap- 
proach their corresponding Markovian counterparts more 
rapidly as the value of the temperature becomes higher. 
Similar to the case of increasing temperature, the time- 
dependent non-Markovian coefficients with large cut-off 
energies saturate to their corresponding Markovian coun- 
terparts more quickly as the bias voltage is increased as 
shown in Figs. |4] and El In summary, the characteris- 
tic times for the non-Markovian behaviors of the time- 
dependent coefficients fpi B \ to differ from their Marko- 
vian counterparts are usually longer for smaller cut-off 
energies, bias voltages and temperatures. 

Next we discuss the time-dependent coefficients of the 
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FIG. 7. (Color online) Time-dependent coefficients of the un- 
conditional master equation: (a) decoherence coefficient D e (t) 
in unit of Po'J^ 1 , (b) damping coefficient y e (i) in units of 7e' / , 
(c) diffusion coefficients h e {t) in units of fvy™ and (d) fre- 
quency renormalization shift [u> 2 (t) + u 2 ] in units of uJoJe' 
for different values of the finite Lorentzian cut-off energy: 
(Ae/ftu ) = 1 (in blue solid lines), 10 (in green dashed lines), 
and 100 (in red dot-dashed lines). The Markovian WBL cases 
are plotted in gray dotted lines. The time is in units of w„ ■ 
The inset in (a) shows the large- value behavior of D e (t) in the 
short-time regime for the case of A e = 100huj o , and the inset 
in (d) shows the long-time behavior of the frequency renor- 
malization shift [u>e(t) + Wj] for the case of A e = hu . Other 
parameters used are eV = 0.1ftu o , and ksT — 0.1huu o . 

unconditional master equation defined in Eqs. (|38 p -(|4i p . 
We plot in Figs. 0(a), (b), (c) and (d) the time-dependent 
decoherence coefficient D e (t), damping coefficient 7 e (t), 
diffusion coefficients h e (t) and frequency renormalization 
shift [o)g (t) + u>c] respectively, for different values of the 
Lorentzian cut-off energy A e . The time-dependent coef- 
ficients are affected primarily by the values of cut-off en- 
ergy A e for small fixed values of bias voltage eV — 0.1fruj o 
and temperature ksT = 0.1hui o , and approach to their 
long-time-limit values in a time scale of about H/A e . 
The time-dependent coefficients with large cut-off en- 
ergies saturate at their corresponding Markovian WBL 
values (in dotted lines), while the coefficient with a small 
cut-off energy of A e = Hlo approach a long-time value 
different from the Markovian WBL value. The contribu- 
tion to the frequency normalization uj 2 (t) comes from the 
imaginary part of the combination of the time-dependent 
tunneling coefficients [see Eq. (|38l) ] or the combination 
of the QPC reservoir correlation kernels (functions). We 
can see from Fig. Eld) that the frequency renormaliza- 
tion shift [£>g(t) + uo 2 ] approaches zero at large times 
since the counterterm frequency contribution us 2 com- 
pensates the frequency renormalization at large times, 
i.e., lo 2 = — (uj^ 1 ) 2 . The diffusion coefficient h e (t) com- 
ing from the contributions of the imaginary part of the 
combination of the tunneling coefficients [see Eq. (|41[) ] 
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FIG. 8. (Color online) Time evolutions of the dynamical vari- 
ables of the NMR for different values of the finite Lorentzian 
cut-off energy: (A e /hu ) — 1 (in blue solid lines), 10 (in green 
dot-dashed lines), and 100 (in red dashed lines), (a) position 
{x(t)) in units of xo, (b) momentum (p(t)) in units of po, (c) 
second moment position (x 2 (t)) in units of Xo, and (d) sym- 
metrized second moment position-momentum {{x,p}(t)) in 
units of xopo- The Markovian WBL cases are plotted in gray 
dotted lines. The NMR is initially in a coherent state with 
{x(0)) = 5xo and (p(0)) = 0. The time is in units of uJo 1 - 
The inset in (c) shows the long-time behavior of (c) for small 
values of {x 2 (t)}. Other parameters used are eV = 0.1foj o , 
k B T = 0.1fiw o , and ~if = 0.12cj o . 



or the imaginary part of the combination of the reservoir 
correlation kernels or functions (or coming from the con- 
tributions of the Cauchy principle values in the Marko- 
vian case) is typically very small compared to the other 
coefficients (see Fig. [7]) and thus is often neglected in the 
discussion of the reduced dynamics of the NMR. We will 
show later that a transient current term coming also from 
the contributions of the imaginary parts of the reservoir 
kernels with a different combination has, however, a sub- 
stantial value and should be kept in order to describe 
correctly the measured time-dependent current. 

Figure [5] shows the numerical results of the dimen- 
sionless mean and covariance values of the dynamical 
variables of the NMR for different values of the cut-off 
energy. The dimensionless mean values (x(t))/xo and 
(p(t))/po oscillate with a frequency of about the NMR 
renormalized frequency oj while the variances {x 2 (t))/xQ 
and (x(t)p(t) + p (t) x(t))/xopo oscillate with twice the 
frequency. As expected, the non-Markovian results with 
larger cut-off energies A e are closer to their Markovian 
WBL results, and there are considerable differences be- 
tween the Markovian WBL cases and the non-Markovian 
cases with the cut-off energies A e comparable to the NMR 
frequency u) a . We note here that our results include the 
frequency renormalization in the non-Markovian cases. 
We can see from Fig. [7^d) that the non-Markovian cases 
with larger cut-off energies have a slightly larger renor- 
malized physical frequency in the short-time region. As 
a result, the initial oscillatory behaviors of (x(t))/xo and 
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FIG. 9. (Color online) Individual contributions of the aver- 
age current I(t) of Eq. (|74|) for Markovian WBL case (in red 
dashed lines) and the non-Markovian cases with a small cut- 
off energy of A e = hw (in blue solid lines). The first term 
in Eq. (|74[) . i.e., Eq. (|75jl . is divided into three parts: (a) the 
isolated QPC tunneling current, (b) the current proportional 
to (x(t)) y and (c) the current proportional to (x 2 (t)}. The 
last three terms of the average current, Eq. H74[) . are plotted 
in (d), (e) and (f), respectively. Each individual contribution 
of the average current is in units of I^f Q — e 2 VGo and the 
time is in units of lo^ 1 . The inset in (c) shows the long-time 
behavior of (c) for small values of Ixx(t). Other parame- 
ters used are eV = 0.1hcj o , ksT — O.lhujo, G x xo/Go = 0.04, 
GpXo/Go = 0.01 , G xx xl/G = 0.002, and 7 f = 0.12w o . 



(p(t))/po with larger A e in Figs.[5Ja) and (b) are slightly 
lower (i.e., with a slightly larger frequency) in the short- 
time region than that of the Markovian WBL case (in 
dotted line). Furthermore, one can also observe that 
the dynamical variables with larger cut-off energies A e 
approach to their steady-state values faster than those 
with smaller cut-off energies, and the steady-state values 
of the dynamical variables in Fig.[S]match well the values 
of the analytical expression of Eqs. (frJ5)) (T7U|) and (j72| . 

Figure [5] shows the differences in time evolutions of 
the individual contribution terms of the average current, 
Eq. (174]) , between the Markovian WBL case and the non- 
Markovian cases with a small cut-off energy of A e = hw . 
We further divide the first term in Eq. |(7S|). i.e., Eq. ((751 . 
into three parts: the isolated QPC tunneling current, 
the current proportional to (x(t)), and the current pro- 
portional to (x 2 (t)). These three parts are plotted in 
Figs-Oa), (b) and (c), respectively. The last three terms 
of the total current, Eq. (f74|) . are plotted in Figs. (9]Jd), 
(e) and (f), respectively. There are considerable differ- 
ences between the non-Markovian and Markovian WBL 
results. In particular, the tunneling coefficients in the 
Markovian WBL are constants in time, and thus the co- 
efficients Aj in front of the individual contributions to 
the total average current become also time-independent. 
As a result, the initial values of the individual contri- 



bution terms of the average current will depend only 
on the initial values of the dynamical variables of the 
NMR. For example, we choose an initial state such that 
(x(0)} ^ and (x 2 {0)) ^ 0. Then the Markovian WBL 
current contributions Ix (t) and Ixx if) start at a finite 
value at time t = 0, i.e., the QPC responses instanta- 
neously to the motion of the NMR to generate a finite 
current at t = (see the red dashed lines in Fig. [5]). 
This is of course not physical. In contrast, in the non- 
Markovian case, the tunneling coefficients and Aj(£) are 
time-dependent and their initial values at the moment 
t = when the QPC detector is brought to interact with 
the NMR are zeros. Thus the individual contribution 
terms of the average current start from zeros (see the blue 
solid lines in Fig. [S]) and will approach their Markovian 
(finite-bandwidth) counterparts at a time scale of h/A e . 
The extra transient current term Iix,P}{t) of Eq. |86() , 
plotted in Fig. (H{f), is proportional to the expectation 
value of the symmetrized product of the position and 
momentum operators of the NMR ({x,p}(t)) and thus 
oscillates with twice the frequency of 0J o . This additional 
term, Eq. l(5rj() . that was generally ignored in the studies 
of the same problem in the literature^— has a coeffi- 
cient proportional to the imaginary part of the combi- 
nation of the QPC tunneling coefficients or bath corre- 
lation kernels (functions) of Im[^ — £§], where and 
£2 are defined in Eqs. 1(13]) and |46|). respectively. Re- 
call that the frequency renormalization ui 2 (t) of Eq. ((55)) 
and the diffusion coefficient h e (t) of Eq. ((41"j) are also 
proportional to the imaginary part of the QPC tunnel- 
ing coefficients or bath correlation kernels (functions) but 
with different combinations. The value of the frequency 
normalization Cb 2 (t) oc Im[£"(i) increases as A e 

increases, and diverges as A e — > 00. So a counter term 
is introduced for the purpose of frequency regularization. 
On the other hand, the value of the diffusion coefficient 
h e (t) oc Im [£f (t) — £f (t)) decreases as A e increases even 
though the individual terms of the imaginary parts of 
£f(i) and £f(i) defined in Eqs. (O and gS]) diverge as 
A e — s- 00. The typical values of h e (t) are however very 
small as compared to other time-dependent decoherence 
and dissipation coefficients (see Fig. [7J and thus h e (t) 
is often neglected. Similar to h e (t), the extra transient 
current I {x ,p}(t) oc [A 8 (t) - A 7 (t)] cx Im [£f(t) - 
decreases as A e increases even though the individual 
terms of the imaginary parts of £°(t) and £f (t) defined 
in Eqs. (|4"4l and (|4"r!)) also diverge as A e — > 00. However, 
I{x,p}(t) has a considerable magnitude and should be in- 
cluded into the time-dependent current. Iix,p}(t) is also 
proportional to ({x,p}(t)} which vanishes in the steady 
state, so I{x.p} {t) exits only in the transient regime. 
As mentioned that I{ X ,p}(t) decreases as A e increases, 
the contribution of Iix,p}(t) to the transient current is 
very small for large cut-off energies. Indeed, we can see 
from the red dashed line in Fig.[9jf) that I{x,p}(t) does 
not contribute to the transient current in the Markovian 
WBL case. 

For clarity, we plot in Fig. [TUf a) the difference between 
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FIG. 10. (Color online) (a) Difference between the total cur- 
rent and the isolated QPC tunneling current, Itot(t) — ii so (t), 
and (b)current in (a) with the contribution of I{x,p) if) of 
Eq. Q86p furthermore deducted for different values of the 
Lorentzian cut-off energy: (A e /fiw ) = 1 (in blue solid 
lines), 10 (in red dashed lines), and 100 (in green dot-dashed 
lines). The Markovian WBL cases are plotted in gray dot- 
ted lines. The time-dependent average current is in units of 
-^iso — e 2 VGo and the time is in units of uj^ 1 - Other parame- 
ters used are eV = 0.1Hlj o , ksT — 0.1huj o , G x xo/Go = 0.04, 
Gpxo/Go = 0.01 , G xx x 2 /Go = 0.002, and 7 f = 0.12w o . 



the total current and the isolated QPC tunneling current, 
i.e., I(t) — Ii S o(t), for different values of the cut-off energy. 
For comparison, we further deduct the contribution of 
I{X,p\(t) of Eq. ([56]) from I (t) — I iso (t) , and the resultant 
current is plotted in Fig. [TUT b). In the steady state, the 
expectation values (x)^^ = (p)^^ = ({^p})^^ = 
0. Thus the total steady-state average current that ap- 
proaches its Markovian long-time value becomes 

I(t ->• oo) = I iso {t -> oo) + Ixx{t -> oo) + I QM (t -> oo). 

(90) 

It has been discussed in Ref. |2l| that in the limit of small 
bias voltages and temperatures (i.e., eV,ksT <C Hlj ), 
the quantum correction current Iqm (t — > oo) in the 
steady state has to cancel Ixx(t ~> oo) as the voltages 
and temperatures are too small to excite the NMR. As 
a result, the steady-state average current in this case is 
equal to that of an isolated QPC junction, Ii so (t —> oo). 
This is indeed the case for the low values of voltage 
and temperature chosen in Figs. [9] and [101 We can 
see from Fig. [H^e) and the inset of Fig. E^c) that the 
steady-state ~ > oo) does cancel the steady-state 
Ixx{t — > oo). As a result, the total average current dif- 
ference I(t) — ii S o(i) in Figs. [T0l vanishes in the steady 
state. We can also see from Figs. [TUTa) and (b) that 
for a large cut-off energy or bandwidth of A e = IOOHljq, 
the evolution of the time-dependent current approaches 



that of the Markovian WBL case closely. Without in- 
cluding the contribution of I{x,p] (t) of Eq. (|86p to the 
average current, there are still substantially quantita- 
tive difference between the Markovian WBL current and 
the non-Markovian currents with finite cut-off energies 
as shown in Fig. UOf b). Furthermore, both significantly 
qualitative and quantitative differences between the non- 
Markovian currents and the Markovian WBL current in 
the short-time region are clearly observed in FigflWa). 
The non-Markovian transient currents with small values 
of cut-off energy are characterized by oscillations with 
large amplitudes and twice the NMR renormalized fre- 
quency as compared to the Markovian WBL one. We 
note that although a significant difference between the 
non-Markovian and Markovian WBL currents Ixx(t) 
in Fig. M, c ) can be observed, the values of Ixx [t) are 
small compared to those of the other individual current 
contribution terms. So this difference in Ixx(t) may 
not be easily identified in the total average current as 
that of I{x,P}(t)- This extra contribution of I{x,p}{t) 
of Eq. ([86]) was completely neglected in the discussion 
of the Markovian current for the same problem in the 
literature^—. We find, however, that this extra sig- 
nificant contribution of I{x,p}(t) in the transient current 
may serve as a witness or signature of the non-Markovian 
features for the coupled NMR-QPC system with finite 
cut-off energies (bandwidths). 



B. Inclusion of the effect of thermal bosonic 
environment 

So far, we do not include the effect of the thermal 
bosonic environment into our numerical calculations. A 
natural question is whether including the effect of the 
bosonic environment changes the picture that the tran- 
sient current can be used to witness non-Markovian ef- 
fects of the coupled NMR-QPC system. Notice that most 
of the individual time-dependent current terms depend 
on the product of two factors: the combination of time- 
dependent coefficients A, (t) , and the time-dependent dy- 
namical variables of the NMR. The inclusion of the effect 
of the thermal bosonic environment affects only the time 
evolution of the NMR dynamical variables but not the 
time-dependent coefficients Ai(t). So if the coupling of 
the NMR to the thermal environment is small compared 
with the coupling to the QPC reservoirs, then the main 
non-Markovian feature in the QPC transient current will 
remain. But if the coupling of the NMR to the thermal 
environment is comparable to the coupling to the QPC 
reservoirs, then the dynamical variables will reach their 
steady-state more quickly. Figure [TT] shows the time evo- 
lutions of (x(t)), (p(t)), (x 2 (t)), and ({x,p}(t)) for differ- 
ent values of the ratio of (7^/7^) that characterizes the 
coupling strength of the NMR-thermal bath relative to 
that of the NMR-QPC reservoirs. One can see in Fig. [11] 
that as the value of the ratio (7^/7^) increases for a 
fixed value of 7( F = 0.12w o , the oscillation amplitudes 
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FIG. 11. (Color online) Time evolutions of the dynamical 
variables of the NMR (a) position {x(t)) in units of xo, (b) 
momentum (p(t)) in units of po, (c) second moment position 
(x 2 (t)) in units of Xq, and (d) symmetrized second moment 
position-momentum ({x,p}(t)) in units of xopo for different 
values of the ratio of 7 D /7<T = 0.001 (in blue solid lines), 0.1 
(in green dot-dashed lines), and 1/3 (in red dashed lines), and 
1 (in gray dotted line). The NMR is initially in a coherent 
state with (x(0)) = 5xo and (p(0)) = 0. The time is in units of 
uj' 1 . Other parameters used are eV = 0.1hu o , ksT — 0.1Huj o , 



--^/y^O.001 

'o 'e 

-7 M /Y M =l/3 



Y^/Yf=0.001 

■ — ^=0.1 
- - -7^^=1/3 

'o 'e 




FIG. 12. (Color online) (a) Difference between the total cur- 
rent and the isolated QPC tunneling current, itot(t) — Iiso(t), 
and (b) current in (a) with the contribution of I{x.p} it) of 



A e = Aq = hw , and 7 e 



0.12w o . 



Eq. (1861) furthermore deducted for different values of the ratio 
of hf = 0.001 (in blue solid lines), 0.1 (in green dot- 
dashed lines), and 1/3 (in red dashed lines), and 1 (in gray 
dotted line). The time-dependent average current is in units 
of I?£ = e 2 VGo and the time is in units of cj^ 1 . Other param- 
eters used are eV — 0.1ftw o , ksT = 0.1ftu o , A e = Ao = hw , 
G x x /G = 0.04, G p x /G = 0.01 , G xx x 2 /G = 0.002, and 
7e M = 0.12lj o . 



of (x(t)) 7 (p(t)), (x 2 (t)), and ({x,p}(t)) diminish and the 
time at which the steady state is reached shifts to the 
short time region since the total damping coefficient be- 
come larger. As a result, the differences in oscillation am- 
plitudes between I(t) — Ii so (t) and I(t) — Ii so (t) — I{x,p}(t) 
become small and the time intervals where the differ- 
ences exist with characteristic oscillation frequency of 
2uj also become shorter (see Fig. IT2)) . The oscillations 
will quickly reach their steady-state values and differ- 
ences will become unobservable if the ratio of (7^/7^) 
becomes much larger than one. 

Another question is whether the case where there are 
no non-Markovian effects in the NMR-QPC system, but 
there are non-Markovian effects induced by the bosonic 
environment, will result in similar non-Markovian fea- 
tures in the transient current. The answer to the question 
can be found as follows. In our simple Lorentzian spectral 
density model, no non-Markovian effects in the NMR- 
QPC system implies that the cut-off energy or bandwidth 
of the QPC reservoir spectral density is very large (i.e., 
A e 3> huj or in the WBL) since the Markovian limit can 
be justified by this condition. In the Markovian case, all 
the coefficients of Aj(t) become time-independent. As 
mentioned previously, the time-dependent coefficient of 
[Ag — A 7 ] oc Im[£° — Q] °f the extra transient current 
I{x,p}{t) decreases as A e increases and becomes very 
small in the (Markovian) WBL. So I{x.p} (t) is very small 



in the Markovian WBL [see, e.g., Fig. HJf)] even though 
({x,p}(t)) has a considerable amplitude in the transient 
regime [see, e.g., Fig. HJd)]. Therefore, if there are no 
non-Markovian effects in the NMR-QPC system (i.e., in 
the WBL), the extra transient current term I{ x p y(t) does 
not contribute even though there may be still significant 
oscillation amplitudes in ({x,p}(t)). Figure [T51 shows the 
current difference of I(t) — Ii so (t) for different values of 
the ratio of (7^/7^) in the case where the QPC reser- 
voirs are in the Markovian WBL (i.e., A e —> 00). The 
non-Markovian feature of oscillations with frequency of 
2lj in I(t)—I iso (t) is unobservable in Fig.Q2]as I{x,p} (t) 
does not contribute and Ixx{t) is too small. The slight 
differences in I(t) — Ii so {t) between the case of A Q = 1 
and A„ = 500 are primarily due to the differences in 
(x(t)) and (p(t)) induced by the non-Markovian bosonic 
environment for the two different values of the spectral 
density frequency bandwidth A Q . 



VIII. CONCLUSIONS 

In summary, we have derived second-order time- 
local (time-convolutionless) non-Markovian conditional 
(n-resolved) and unconditional master equations of the 
reduced density matrix of a NMR subject to a measure- 
ment by a low-transparency QPC or tunnel junction de- 
tector and an influence by a thermal environment. Our 
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FIG. 13. (Color online) Current difference of I(t) - 7 iso (t), 
for different values of the ratio of (a) jo 1 /je 1 = 1/3 and 
(b) -jo' he 1 = 0.001. The QPC reservoirs are in the WBL 
(i.e., A e — > oo), while the frequency bandwidth of the spec- 
tral density of the thermal bosonic bath is Ao = uio (in blue 
solid lines) and Ao = 500a>o (in red dashed lines). The time- 
dependent average current is in units of I{f a = e 2 VGo and 
the time is in units of lo^ 1 . Other parameters used are eV = 
0.inw o , k B T = 0.1fiw o , G x x /G = 0.04, G p x /G = 0.01 , 
G xx Xo/G = 0.002, and = 0.12w o . 

non-Markovian master equations implemented with the 
reservoir memory correlation prescription going beyond 
the WBL allow us to study the memory effect of the 
non-equilibrium QPC fermionic reservoir and the equi- 
librium bosonic thermal bath on the NMR. Our non- 
Markovian master equations with time-dependent coeffi- 
cients reduce, in appropriate limits, to various Markovian 
versions of master equations in the literature. Further- 
more, our non-Markovian master equations are valid for 
arbitrary temperatures of the thermal environment and 



QPC reservoirs (detector), and for arbitrary bias voltages 
as long as the perturbation theory up to the second order 
in the system-detector and system-environment coupling 
strength holds. 

We have found considerable differences in dynam- 
ics between the non-Markovian cases and their Marko- 
vian counterparts. The fact that the QPC detector in- 
duces a back action on the NMR and the motion of the 
NMR modulates the current through the QPC are taken 
into account self-consistently. We have also calculated 
the time-dependent transport current through the QPC 
which contains information about the measured NMR 
system. We have found an extra transient current term 
°f I{x.p}(t) °f Eq. ([86]). This extra term, with a co- 
efficient coming from the combination of the imaginary 
parts of the QPC reservoir correlation functions, was gen- 
erally ignored in the studies of the same problem in the 
literature. But we find that it has a substantial contribu- 
tion to the total transient current in the non-Markovian 
finite-bandwidth case and differs qualitatively and quan- 
titatively from its Markovian WBL counterpart. Thus 
it may serve as a witness or signature of non-Markovian 
features for the coupled NMR-QPC system. 
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